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to . Abstract 

For the GMRT, we calculate the expected signal from redshifted HI emission 
^ ■ at two frequency bands centered at 610 and 325 MHz. The study focuses 

on the visibility-visibility cross-correlations, proposed earlier as the optimal 
statistical estimator for detecting and analyzing this signal. These correla- 
tions directly probe the power spectrum of density fluctuations at the redshift 
where the radiation originated, and thereby provide a method for studying 
the large scale structures at large redshifts. We present detailed estimates 
of the correlations expected between the visibilities measured at different 
baselines and frequencies. Analytic fitting formulas representing the salient 
features of the expected signal are also provided. These will be useful in 
planning observations and deciding an optimal strategy for detecting this 
signal. 
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1 Introduction 



In two earlier papers (Bharadwaj, Nath & Sethi, 2001; hereafter referred 
to as Paper I and Bharadwaj & Sethi, 2001; hereafter referred to as Paper 
II) we have studied the possibility of detecting redshifted 21 cm emission 
from neutral hydrogen (HI) at high redshift sand using these observations to 
probe large scale structures. Our approach has been based on the fact that 
the HI density in the redshift range 1 < z < 3.5 is known from observations 
of Lyman-a absorption lines seen in quasar spectra. These observations 
currently indicate Q gas (z), the comoving density of neutral gas expressed as 
a fraction of the present critical density, to be nearly constant at a value 
Qgas{ z ) ~ 10~ 3 for z > 1 (Peroux et al. 2001). The bulk of the neutral gas 
is in clouds which have HI column densities greater than 2 x 10 20 atoms/cm 2 
(Peroux et al. 2001, Storrie-Lombardi, McMahon, Irwin 1996, Lanzetta, 
Wolfe, & Turnshek 1995). These high column density clouds are responsible 
for the damped Lyman-a absorption lines observed along lines of sight to 
quasars. The flux of HI emission from individual clouds (< lOyuJy) is too 
weak to be detected by existing radio telescopes unless the image of the 
cloud is significantly magnified by an intervening cluster gravitational lens 
(Saini, Bharadwaj and Sethi, 2001). Although we may not be able to detect 
individual clouds, the redshifted HI emission from the distribution of clouds 
will appear as a background radiation in low frequency radio observations. In 
the two earlier papers, and in the present paper we investigate issues related 
to calculating the expected signal and detecting it. 

The hope of detecting this signal lies in the fact that the distribution of 
the HI clouds is expected to be clustered. The clustering pattern of the HI 
clouds will be reflected in the angular and spectral distribution of the the 
redshifted 21 cm emission. It is the aim of this paper, and the two previous 
papers, to characterize the expected properties of this signal and to pro- 
pose a optimum statistical estimator to extract this signal from the various 
foregrounds and the system noise which will also be present in any low fre- 
quency observation. In paper II we proposed the cross-correlations between 
the visibilities measured at different baselines and frequencies as a statistical 
estimator to optimally detect the HI signal. This has the advantage that 
the system noise contribution to visibilities measured at different baselines 
and frequencies are uncorrelated. Further, in Paper II we showed that the 
visibility-visibility cross-correlations directly probes the power spectrum of 
density fluctuations at the epoch where the HI emission originated. In this 
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paper we present results of a detailed investigation of the signal expected for 
this statistical estimator for observations with the GMRT. 

We next present a brief outline of this paper. In section 2 we review 
the formalism used to calculate the visibility-visibility cross-correlations. In 
section 3. we present the values of a set of useful parameters to do with 
the GMRT (Swarup et at. 1990), the cosmological model chosen for these 
investigations and the HI distribution. In Section 4. we present the results 
of our investigations, and in Section 5 we discuss the results and present 
conclusions. 

2 Formalism. 

We consider a situation where an interferometric array of N antennas is used 
to carry out low frequency radio observations. The observations are carried 
out at NC frequency channels {ui, vi, is 3 , v^c} covering a frequency band 
B centered at the frequency v c . For the purpose of this paper we assume that 
the antennas are distributed on a plane, and the position of each antenna is 
denoted by a two dimensional vector dj. The antennas all point in the same 
direction along the unit vector n which we take to be vertically up wards. The 
beam pattern A(9) quantifies how the individual antenna, pointing upwards, 
responds to signals from different directions in the sky. This is assumed to be 
a Gaussian A(9) = e~ e with 9 <C 1 i.e. the beam width of the antennas 
is small, and the part of the sky which contributes to the signal can be well 
approximated by a plane. 

The quantity measured in interferometric observations is the visibility 
V(U, v) which is recorded for every independent pair of antennas (baseline) at 
every frequency channel in the band. For any pair of antennas, the visibility 
depends on the vector d = dj — dj joining the position of the two antennas. It 
is convenient to express the visibility as a function of the variable U which is 
d expressed in units of the wavelength i.e. U = d/A. Typically, the visibility 
is recorded for U in a range ?7 min <| U |< ?7 max . The visibility corresponding 
to U = is the flux measured by a single antenna, and this is usually not 
recorded in an interferometric array. 

In paper II we calculated the HI contribution to the visibility V(U, v). 
This can be expressed as 



V{U, u)=I, 
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The frequency v appears in many of the terms on the right hand side 
of equation (1) to indicate that these are to be evaluated at the redshift 
z = v e j v — 1 where the radiation originated. The first term l v is the specific 
intensity of the redshifted HI emission expected if the HI were uniformly 
distributed. This is given to be (from eq. (3) of Paper II) 



2.7 Jy 

degree 2 



hQ gSuS (z) 



Hp 
H(z) 



(2) 



where H(z) is the Hubble parameter, and Hq is its present value. 

It should be noted that eq. (14) of Paper II is incorrect. Equation 3 
of paper II is correct. The error came about when expressing nm(z), the 
comoving number density of HI atoms in terms of Qm- As a result of this 
and an error in the computer programs that were used, the predictions for 
610 MHz in Paper II are also incorrect. 

The HI emission originates from the gas located at a distance r v in the 
direction n. The term A^j(k, v) in equation (1) is the Fourier transform of 
the density contrast of the HI distribution in redshift space. Any Fourier 
mode k can be decomposed into two parts, a component k\\ parallel to n 
and a component perpendicular n. Assuming that the HI density fluc- 
tuations evolve according to linear theory, and that these may be related to 
the fluctuations in the underlying dark matter distribution through a linear 
bias parameter b(z), we have 



A s m (k,u) =b(z)D(z) 



1 + /?(,)§ 



A(k) 



(3) 



where D(z) is the growing mode of density perturbations (Peebles, 1980), 
P( z ) = f(Qm,QA)/b( z ) (Lahav et al. 1991) is the linear redshift distortion 
parameter, and A(k) is the Fourier transform of the density fluctuations of 
the dark matter distribution at the present epoch. 
We also have 



(A(k)A (k )) = (27r) 3 ^(k - k )P(k) 



(4) 



where () denotes ensemble average, and P(k) is the power-spectrum of dark 
matter fluctuations at the present epoch. 

The quantity a u (XJ) in equation (1) is the Fourier transform of A(9), the 
beam pattern. For a Gaussian beam, the Fourier transform also is a Gaussian 



4 



and we have 



ixQl v exp 
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This function is sharply peaked around = 27rU/r„, and only values of k^ 
in the interval | Akj_ |< l/(r v 6Q V ) centered around this value will contribute 
to the integral in equation (1). It follows that the visibility V(U, v) picks up 
contributions only from those fluctuations A| II (k, v) for which kj_ 2ir\J/r u . 

A point to be noted is that in the Gaussian model for the beam pattern 
its Fourier transform is also a Gaussian which is a non-compact function. In 
reality, the Fourier transform of the beam pattern is compact and it is zero 
outside a radius D/X i.e. a v (U - ^) = for | U - ^ |> D/X. This 



fact plays a role in the later discussion and we shall refer back to it at the 
appropriate place. 

In this paper we shall study in some detail the cross-correlation (V(Ui, za)U*(U 2 , ^2)) 
between the visibilities measured at two different baselines and frequencies. 
Using equations (1), (3), (4), this turns out to be 



(y(Ui,i/i)y*(u 2 ,^)> 
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the imaginary part being zero. We shall use equation (6) to calculate the 
correlations between the visibilities at different baselines and frequencies. We 
next discuss a few approximations which can be used to simplify equation (6) 
and which help in interpreting this equation. We test the range of validity 
of these approximations and comment on this later in the paper. 

Typically, the bandwidth B is small compared to the central frequency 
v c , and the quantities I u , b(z), D(z), f3(z), 6 0u and r u in equation (6) do 
not vary substantially across the band. We can evaluate these at v c instead 
of calculating them separately at v\ and u 2 , and at some instances we do 
not explicitly show the subscript v\ or v 2 for these quantities. The term 
cos[/c||(ri — r 2 )] is an exception because — r 2 ) may be large even if 

Au = u2 — v 1 <^ u c . This term may be approximated as cosf/syr^Ai/], where 
r v = dr u / du. 
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The term ai(Ui — ^^(L^ — ^r) in equation (6) is a product of two 
functions, one sharply peaked at kj_ = 27rUi/r and another at kj_ = 27tU2/t. 
The two peaks have very little overlap if | Ui— U2 1/0q, and the visibilities 
measured at two such baselines will be uncorrelated. There will be correla- 
tions only between visibilities measured at different frequencies for nearby 
baselines, ie. | Ui — U 2 |< l/9 . For the situation where Ui = U 2 3> l/#o, 
we may use the approximation (Paper II) 



'2tt 3 9 2 
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(7) 



whereby equation (6) is considerably simplified and we have 



(V(U, v)V*(U,v + Av)) 
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where k = ^J(2%U/r) 2 + fcjj. 

Equation (8) provides considerable insight on the visibility- visibility cross- 
correlation in the situation where Ui = U 2 = U. The point to note that the 
visibility signal measured at a baseline U receives contributions only from the 
fluctuations A(k) for which kj_ = (2 7rU)/r, ie. h± is fixed but k\\ can have 
any value. As a consequence the correlation (V(U, z/)l / *(U, u + Au)) directly 
probes the power spectrum P(k) at all Fourier modes k > k min = 2nU/r. 

In the next sections of this paper we use the formulas presented here 
to investigate the HI signal expected for the GMRT and study how this is 
related to the large-scale structures at the redshift where the HI emission 
originated. We also compare the predictions of eq. (6) with those of (8) and 
test the range of validity of the approximations discussed above. 



3 Some Useful Parameters 
3.1 GMRT 

The Giant Meterwave Radio Telescope (GMRT) has 30 radio antennas of 45m 
diameter each. We present below, in tabular form, some of the parameters 
relevant for the proposed HI observations. We have restricted our analysis 
to only two of the GMRT frequencies which correspond to HI in the range 
1 < z < 3.5. 



6 



Frequencies v c (MHz) 


610 325 


Possible Bandwidths - B (MHz) 


16 8 4 2 1 0.5 


No of Channels - NC 


128 


Antenna Separations - d (m) 


d mhl rsj 60 d max ~ 25 x 10 3 





^min 


^max 


8q ~ 0.6 X $FWHM 


610 


123 


51 x 10 3 


0.54° 


325 


65 


27 x 10 3 


1.08° 



3.2 Cosmological Model 

For the background cosmological model we have used H Q = 100/ikm/s/Mpc 
with h = 0.7, Q m0 = 0.3, Q\ = 0.7, n Bary ono = 0.015 h~ 2 and Q gas = 
1. x 10~ 3 . The power spectrum of dark matter density fluctuations is nor- 
malized to COBE (Bunn & White 1996), and its shape is determined using 
the analytic fitting form for the CDM power spectrum given by Efstathiou, 
Bond and White (1992). The value of the shape parameter turns out to 
be T = 0.2 for the set of cosmological parameters used here. The power 
spectrum is shown in Figure 1. 



100000 
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Figure 1: This shows the part of the power spectrum probed by the smallest 
baselines in GMRT for the two frequencies indicated in the figure. The power 
spectrum is normalised to COBE and is shown for z = 0. 
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Very little is known about the bias b which relates the HI fluctuations to 
the dark matter fluctuations, and we have used a fixed value 6=1 through- 
out. 

3.3 HI. 

We first present the values of the redshift, comoving distance and the ex- 
pected specific intensity of the HI distribution for the two GMRT frequencies 
of our interest. We also show values of r = dr/du which are useful in calcu- 
lating the comoving separation corresponding to a given frequency interval 
ie. Ar = r Au. The last column in the table below shows the value of /c min 
for the baseline U = 100. This represents the smallest Fourier mode probed 
by this baseline. 





z 


r (Mpc) 


/ (Jy/deg 2 ) 


r (Mpc/MHz) 


k min (Mpc" 1 ) 


610 


1.33 


4030 


9.0 x 10~ 4 


7.7 


0.16 


325 


3.37 


6686 


3.7 x 10~ 4 


11.3 


0.09 



In Figure 1 we show the part of the power spectrum which will be probed 
by U m [ n the smallest baseline available on GMRT. The value of k m - m will be 
higher for the larger baselines and these will probe a smaller part of the power 
spectrum. 



4 Results 

4.1 Single Dish Observations 

It is possible to use the GMRT as a set of 30 single dish antennas instead of an 
interferometric array. In such observations the flux incident on each antenna 
is recorded and the signal from the 30 antennas is combined. The flux F{y) 
is the visibility V(XJ, u) at U = 0, and the flux expected from the redshifted 
HI emission can be calculated using equation (1). We use equation (6) to 
calculate the correlation between the flux measured at different frequencies 
(F(u)F(u + Au)). The results for the two GMRT bands of interest, are 
shown in figures 2 and 3. The approximate formula (eq. 8) is valid only 
when | U |>> l/#0) an d this cannot be used here. 
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1 2 3 4 5 

Av MHz 

Figure 2: This shows the contribution from redshifted HI emission to the cor- 
relation between the flux expected at two different frequencies in the 610MHz 
band for single dish observations using GMRT . The analytic fit (eq. 9) is 
also shown. 

We find that the HI flux at neighbouring frequencies are correlated. The 
correlation signal in the 610MHz band is nearly an order of magnitude higher 
than the correlation signal in the 325MHz band. In both the bands the 
correlation falls exponentially as Av is increased. For Av < 2MHz, the 
correlation can be very well approximated by the fitting formula 

ii?t \w aw n ( ~ Az/ \ f Ceio = 5 x i(T 8 Jy 2 

{F( V )F [V + A,)) = Cexp ( njsi -) | Ca25 = 5 x lQ _ 9 / y2 (9) 

The fitting formula (eq. 9) fails at larger Av where (F(v)F(v+Av)) becomes 
negative indicating that the fluxes have a weak anti-correlation for Av > 
4MHz. 

4.2 Interferometric Observations 

We next consider interferometric observations with GMRT. The correlation 
(V(Ui, vi)V*(XJ2, v 2 )) is maximum when Ui = U2, and we consider this 
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Figure 3: This shows the contribution from redshifted HI emission to the cor- 
relation between the flux expected at two different frequencies in the 325MHz 
band for single dish observations using GMRT . The analytic fit (eq. 9) is 
also shown. 

situation first. In addition, it is a good approximation to represent the cor- 
relation as a function of Av = v 2 — v\, provided v\ and v 2 are in the same 
GMRT band, Our results are shown in figures 4. and 5. 

We have calculated the correlations using two different formulas i. e. equa- 
tions (6) and equation (8). The latter can only be used when Ui = U 2 = U, 
and it is based on an approximation which is valid only when U ^> 1/9q. 
We find that the approximate formula (eq. 8) matches equation (6) quite 
well (~ 10%) at the smallest baseline U min for both 610MHz and 325MHz 
(figures 4 and 5). The agreement is better than 10% for the larger baselines. 

The correlation signal is the strongest at the smallest baseline C/ m i n and 
falls quickly for the larger values of U. The larger baselines receive contribu- 
tions from a smaller part of the power spectrum (figure 1.) and the fall-off 
with U reflects the shape of the power spectrum at the mode k ~ (2irU/r). 
At a fixed U, the visibilities at adjacent frequencies are correlated, the cor- 
relation falls approximately exponentially with increasing Av. At the larger 
baselines the correlation decreases faster with increasing Av. At all baselines, 
the correlation crosses zero and there are anti-correlations for large Av. We 
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Figure 4: This shows the correlation between the visibilities V(U, v) and 
V(U, v + Av) expected for the same baseline U at two different frequencies. 
The results are shown for different values of U (shown in the figure) starting 
from C/ m i n for the 610MHz band. The results have been calculated using 
equation (6), the predictions of the approximate formula (eq. (8) are also 
shown for U mm . 



find that in the range of U and Av where the visibility-visibility correlation 
signal is large, it can be approximated by the fitting formula 

(y(U, v)V*(U, v + Ai/)) = C „exp (=^f-) (10) 

C 610 = 1.6 x 10" 8 Jy 2 M- L2 C 325 = 1.8 x 10 " 9 Jy 2 M~ a9 

where where u = (C//100) and Av is in MHz. 

We next calculate the correlation (V(U, v)V*{\J + A\J , v + Av)) between 
the visibilities at two different baselines and frequencies. We use equation (6) 
to calculate this. In general, the behaviour of the correlation will depend on 
both, the magnitude and the direction of AU. We have separately considered 
two situations (1.) AU parallel to U and (2.) AU perpendicular to U. The 
results are presented in figures 6 and 7 for U = 200. The results are similar 
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Figure 5: This shows the same thing as figure 4. for the 325MHz band. 



for other values of U. 

We find that the correlation falls rapidly as AU increases. The AU de- 
pendence is found to be anisotropic for the 610 MHz band, . The correlation 
falls faster when AU is parallel to U than when it is perpendicular to U. 
The AU dependence is found to be isotropic for the 325 MHz band, In both 
the bands the AU dependence can be well approximated by a Gaussian 



(V(U,v)V*(U + AU,v + Ais)) = exp 



(AU, 



AC/ 1 



&6io — 50, &6io — 53 



\ a u ) \ b u 
(V(U,u)V*(U,u + Au)) 

0-325 = &325 = 28 



X 



(11) 



where AU\\ is the component of AU parallel to U, and AU± is the compo- 
nent perpendicular to U. This fitting function is found to be a very good 
approximation over a large range of of U, AU and Au. 

It should be noted that contrary to the predictions of our fitting formula 
(eq. 11), the cross-correlations between the different baselines is exactly zero 
i.e. (V(U, v)V*(U+AU, u+Au)) = if j AU |> D/X. This is a consequence 
of the fact that we have modeled the Fourier transform of the beam pattern as 
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Figure 6: This shows the expected correlation between the visibilities V(XJ, v) 
and V^U+AU, u+Au) at two different baselines and frequencies. The results 
are shown at U — 200 for different values of AU (shown in the figure), when 
AU is parallel to U (solid line) and AU is perpendicular to U (dashed line). 
These results are for the 610MHz band. 



a Gaussian whereas it actually has compact support (discussed in Section2.), 
The cross-correlations are zero for | AU | > 90 at 610 MHz and | AU | > 50 
at 325 MHz. This discrepancy does not significantly change our results as the 
Gaussin function in equation (11) is extremely small (< 1%) at the values of 
I AU I where the cross-correlations vanish. 



5 Discussion and Conclusions. 

We have investigated in detail the redshifted HI emission signal expected 
in GMRT observations in the 325MHz and the 610MHz frequency bands. 
The properties of the visibility-visibility cross-correlation (V(U, v)V*(XJ + 
AU, v + Av)) which we propose as a statistical estimator for detecting and 
analyzing the signal have been studied. The signal 

we find that (V(U min , v)V*(U min , v)) is ~ 1.3 x 10~ 8 Jy 2 at 610 MHz and 
~ 2.5 x 10~ 9 Jy 2 at 325 MHz. The signal falls at larger baselines, and it drops 
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Figure 7: This is same as figure 6 except that it is for the 325MHz band. 
Here the results do not depend on whether AU is parallel or perpendicular 
to U. 

by an order of magnitude by U ~ 10 x U min . Fourteen of the thirty GMRT 
antennas are located within 1km. These baselines could be used for detecting 
the signal whereas the larger baselines where there is very little signal could 
be used as a control to test that what we detect is actually the HI signal and 
not an artifact. 

The correlation {V(U, v)V*(XJ, v + Av)) fall exponentially with increas- 
ing Av, the width being around 0.5 MHz at U m \ n . The correlations fall faster 
at the larger baselines. The Av dependence of the correlation will be cru- 
cial in detecting the HI signal, for neither the system noise nor the galac- 
tic/extragalactic continuum radiation are expected to have such a feature. 

The correlation (V(U, z/)l / *(U+AU, v+Av)) falls rapidly with increasing 
AU, the AU dependence being well described by a Gaussian. The correlation 
falls by 50% when | AU |~ 42 at 610MHz and | AU |~ 23 at 325MHz. The 
correlation falls by 10% when | AU |~ 16 at 610MHz and | AU |~ 9 at 
325MHz. The AU dependence plays an important role in three different 
contexts in this discussion. 

a. For a pair of antennas at a separation d, the value U = d/A changes as 
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v varies across the frequency band. The correlation between the visibil- 
ities measured by this pair of antennas at two different frequencies will 
be (VTU,i/)V*(U + AU,i/ + Ai/)), where AU = -(Ai//i/)U arises due 
to the change in wavelength. . For U < 500 and Av < 1MHz, which 
is the range where there is significant correlation, we find | AU |< 2. 
Our investigations show that this effect will cause the correlation to 
fall by a factor less than 10% from (V(U, u)V*(U, v + Av)), and we 
can ignore the change in U due to the change in frequency. 

b. The rotation of the earth causes all the baselines to change in time 
and for a time interval At we have AU = AtcJxU. This has to be 
taken into account when deciding on the time interval At over which 
the visibility can be integrated without its value changing significantly. 
Our analysis shows that the correlation in the visibility falls by 10% 
or less if AU < 10 and we use this as a the criteria to determine At. 
In the range U < 500 where we expect most of the signal we obtain 
At ~ 5 min. 

c. As the earth rotates, GMRT observations will produce visibility data 
for a very large number of baselines U. The values of U will not 
be uniformly distributed in the plane of the array, and there will be 
regions which are more densely sampled than others. It will be useful 
to reduce the volume of data by combining the data in the densely 
sampled regions and thereby produce a set of values of the visibility on 
a regular grid of U values. Our analysis shows that the grid has to be 
at an intervalAt/ ~ 40 at 610MHz and ~ 25 at 325MHz, the visibility 
signal gets uncorrelated after this interval. 

Two issues, the system noise and the contribution of the galactic and 
extragalactic continuum sources have not been addressed here. Further, we 
have modeled the HI gas as having a continuous distribution, whereas in 
reality the HI gas is in discrete, small clouds. Investigations are currently 
underway on these issues and the results will be communicated in forthcom- 
ing publications. 

Finally, the visibility-visibility correlation signal depends the equation 
of state of the dark matter and the dark energy at the the epoch when 
the HI emission originated. This dependence comes in through the redshift 
space distortion parameter, the comoving distance and its derivative, all of 
which depend on the equation of state of the universe. We shall address 
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the possibility of using HI observations to probe the equation of state of the 
universe at at high redshifts in a future paper. 

Acknowledgment. SB would like to than Jayaram N Chengalur and Shiv 
K Sethi for useful discussions. SB would also like to acknowledge financial 
support from BRNS, DAE, Govt, of India, for financial support through 
sanction No. 2002/37/25/BRNS. SKP would like to acknowledge the As- 
sociate Program, IUCAA for supporting his visit to IIT,Kgp and CTS, IIT 
Kgp for the use of its facilities. 

References 

Bharadwaj S., Nath B. & Sethi S.K. 2001, JAA. 22, 21 
Bharadwaj, S. & Sethi, S. K. 2001, JAA, 22, 293 
Bunn E. F. & White M. 1996, ApJ, 460, 1071 

Efstathiou, G., Bond, J. R. & White, S. D. M. 1992, MNRAS, 250, lp 

Lahav O., Lilje P. B., Primack J. R. and Rees M., 1991, MNRAS, 251, 
128 

Lanzetta, K. M., Wolfe, A. M., Turnshek, D. A. 1995, ApJ, 430, 435 

Peebles, P. J. E. 1980, The Large-Scale Structure of the Universe , 
Princeton, Princeton University Press 

Peroux, C, McMahon, R. G., Storrie-Lombardi, L. J. & Irwin, M .J. 
2001, astro-ph/0107045 

Saini, T., Bharadwaj, S. & Sethi, K. S. 2001, ApJ, 557, 421 

Storrie-Lombardi, L.J., McMahon, R.G., Irwin, M.J. 1996, MNRAS, 
283, L79 

Swamp, G., Ananthakrishan, S., Kapahi, V. K., Rao, A. P., Subrah- 
manya, C. R., & Kulkarni, V. K. 1990, Curr. Sci., 60, 95 



16 



